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RESUMEN 


En este artículo se presenta un nuevo método de diseño de Modelo de Network Control System (NCS) 
con retardos en la transmisión. Este método parte de un modelado del NCS usando variables de estado en 
el dominio del tiempo, en el que se ha revisado todas las posibles interacciones que pueden aparecer en el 
conjunto operador — maestro— medio— esclavo. El método de diseño se basa en la obtención y resolución 
de un conjunto de ecuaciones que consiguen la convergencia del estado del esclavo y del maestro, 
logrando de este modo que el esclavo siga al maestro, así como establecer la dinámica de dicha 
convergencia y del esclavo. Entre los aportes principales realizados en este artículo se encuentra el 
modelado de un NCS con retardos en variables de estado contemplando todas las posibles interacciones 
que pueden existir entre el operador, el maestro, el esclavo y el entorno. Utilización de la aproximación 
de Taylor de primer orden para obtener un modelo de estado del NCS sin la presencia de términos 
retardados en el tiempo. El método de diseño que logra la convergencia de estados entre el maestro y el 
esclavo a pesar de la existencia de retardos, y además es capaz de fijar tanto la dinámica de la 
convergencia, como la del esclavo, además ell método de diseño se reduce simplemente a la resolución de 
un sistema de ecuaciones para obtener los parámetros de regulación del esquema de control. 


Palabras Clave: Control, Modelado, Network Control System, Simulación, Taylor. 


ABSTRACT 


This paper presents a new design method Model Network Control System (NCS) with transmission 
delays is presented. This method of modeling the NCS using state variables in the time domain, which 
has reviewed all possible interactions that may occur in the joint operator - average slave master. The 
design method is based on obtaining and solving a set of equations that achieve convergence the state of 
the slave and the master, thus achieving the slave to follow the master, and to establish the dynamics of 
this convergence and slave. Among the major contributions made in this paper is the modeling of a NCS 
with delays in state variables considering all possible interactions that may exist between the operator, the 
master, the slave and the environment. Using the Taylor approximation of first order to obtain a state 
model without the presence of NCS delayed terms in time. The design method that achieves convergence 
of states between the master and slave despite the existence of delays, and it is capable of binding both 
the dynamics of convergence, as the slave also ell design method reduces simply to solving a system of 
equations to obtain the control parameters of the control scheme. 


Keywords: Control, Modeling, Network Control System, Simulation, Taylor. 
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L INTRODUCCIÓN 


G. Canata et al. en 1996 propusieron un método de diseño por aproximación de estados, 
método efectivo a lo largo de los años, pero en este artículo se ha utilizado la 
aproximación de Taylor para modelar el retardo existente en la transmisión, lo cual 
constituye una idea original en el modelado y diseño de de Network Control System 
(NCS) con retardos. 


En publicación de Rasmussen en 2005, se describe un problema en el diseño cada vez 
mayor y los sistemas intensivos en software de construcción para la exploración del 
espacio que los métodos tradicionales de ingeniería no se mantienen al día con la 
complejidad de los tipos de sistemas que se están concibiendo. Aunque el diseño 
riguroso se aplica generalmente a la ingeniería de sistemas de hardware y subsistemas, 
lo mismo no se hace generalmente para el software, a pesar de que el software se 
encarga de una parte creciente de los esfuerzos de ingeniería y costo, y el software se 
encarga de la realización de funciones cada vez más importantes en formas cada vez 
más sutiles. 


Su solución a esta situación es una metodología rigurosa ingeniería de sistemas, 
Análisis, y una arquitectura de software orientada al control basado en el Sistema de la 
Misión de Datos de Jet Propulsion Laboratory según su manual estándar (NASA,1998) . 
Este documento presenta relación con esta tesis debido a que se busca una solución a la 
presente situación proponiendo al menos una especie de metodología de ingeniería, 
Análisis, y proponer una arquitectura orientada al control NCS. 


La importancia y novedad del método de diseño propuesto en este artículo radica en la 
posibilidad de la obtención de un esquema de control para que la dinámica del esclavo, 
así como la dinámica de la convergencia de estados del maestro y del esclavo, tengan el 
comportamiento deseado algo diferente a lo expuesto en 1998 por Fujimoto y Tanaka, 
los cuales realizan la convergencia en sistemas bilaterales a través de compensación de 
parada y ganancia (H.S. et al., 1998), a partir del retardo existente en la comunicación, 
de la constante de reflexión de esfuerzos y de las características dinámicas del maestro, 
del esclavo y del entorno. 


En primer lugar se va exponer el modelado en representación interna del NCS con 
retardos en la transmisión. Para explicar el modelado y el método de diseño se 
considerará un NCS de 1 GDL en el que el maestro y el esclavo están representados por 
ecuaciones diferenciales de segundo orden. Se hará especial énfasis en el modelado del 
entorno y en la utilización de la aproximación de Taylor para modelar el retardo. Una 
vez presentado el modelado, se presentará el método de diseño del sistema de control 
mediante convergencia de estados. 


El artículo está estructurado de la siguiente manera: en la sección II se presenta el 
desarrollo de la investigación: variables de estado y polinomio de Taylor Sistemas 
NCS. En la sección MI se presentan los resultados de los experimentos usando el 
software y discusión en relación al experimento, en la sección IV se presentan las 
conclusiones y en la sección V las Referencias Bibliográficas 
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II. DESARROLLO. 


El modelado del NCS maestro — esclavo con retardos en la transmisión se ha realizado 
en representación interna considerando todas las posibles interacciones que pueden 
aparecer en el conjunto operador — maestro — esclavo — entorno (Ambike, 2009). Este 
modelado se muestra en la Figura 1. 


Figura 1. Modelado de NCS con retardos en la 
transmisión 


El sistema maestro y esclavo, respectivamente, se representan en el espacio de estados 
de la siguiente forma. 


Xm(t) = Am*m(t) + Bmum(t) 
Ym(t) = Cm Xm(t) 
Xs(t) = Asxs(t) + Bsus(t) 
ys(t) = cyxs(t) 

dónde: 

* Xm y xs son las variables de estado, 

* Ym Y y, representan las salidas, 

* Am y As son las matrices de estado, 

Bm y Bs son las matrices de entrada, 

y Cm y Cs son las matrices de salida. 

Las matrices que aparecen en el modelo son las siguientes: 

G2: Define la influencia de la fuerza ejercida por el operador sobre el maestro en el 
esclavo. 

Rm: Define la interacción esclavo — maestro. 

Rs: Define la interacción maestro — esclavo. 

Km: Matriz de realimentación de estado del maestro. 

Ks: Matriz de realimentación de estado del esclavo. 


Cabe resaltar que la topología de NCS que se muestra en Figura 1 en la que se 
consideran todas las posibles interacciones que pueden aparecer en el conjunto operador 
— maestro — esclavo — entorno, supone una idea innovadora en el modelado de sistemas 


35 


Ciencia e Ingeniería 2015 Vol. (2) N* (2) Velásquez S, Velásquez R. y Romero H. 


MODELADO DE UN NETWORK CONTROL SYSTEM CON RETARDOS EN LA TRANSMISIÓN 


de control NCS con retardos. 


Análisis del modelo 
Se va a realizar el análisis suponiendo que el maestro y el esclavo constan de un único 
grado de libertad, ver Figura 2. El modelo dinámico de un elemento de un único grado 
de libertad es el siguiente: 


JÓ(O) + b6(t) + mgl sinó (t) = ult) 
dónde: 
J= ml es la inercia del elemento, b es el coeficiente de rozamiento viscoso, l es la 
longitud del elemento, g es la aceleración de la gravedad, m es la masa del elemento, 


0 (t) representa la posición (ángulo girado) en el instante t y u(t) es el par de control 
aplicado en el instante t (H.P., 1999). 


Figura 2. Modelo simplificado de un robot de 1 GDL. 


Para efectuar el análisis, se va a considerar el modelo dinámico simplificado lineal: 
J6(t) + b6(t) = u(t) 


Se va a expresar el modelo dinámico en el espacio de estados mediante las ecuaciones: 
x(t) = Ax(t) + Buít) 
y(t) = Cx(t) 


dónde el vector de estados es de la forma: 


x (6) 
iS Leto 


Para modelar el sistema en el espacio de estados se van a elegir variables de estado con 


sentido físico. En concreto se van a considerar como variables de estado la posición y la 
velocidad del elemento. 


Sea x,(t) =0 (t) y x2 (t) =x,(t) =0(t) , es decir, el estado xl representa la 
posición del elemento, y el estado x2, la velocidad. A partir de la ecuación (6) se va a 


calcular la ecuación de estados (Ferrel, 1986; García et al., 2000 y Lían et al., 2007). 


(6) = x2(6) = 0(t) 
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%2 (6) = X,(6) =Ó(0) 
= Lo) PO 
j j 


ón bs t 
= ya pa 


y(t) = x1(t) 


Por tanto, el modelado en variables de estado con sentido físico del sistema es: 


Leste) xa(£) Luo 


20] =lo - » po, 


yO =b 0 po 


De esta forma, la representación en el espacio de estados del sistema maestro es 
(Coiffet, 1990): 


Anc 


0] =lo - A po, 


Es col Xm2[t) 


Ym(t) =[1 0] eco] * 


La representación en el espacio de estados del sistema esclavo es: 


Lu 


E=0] =lo - a] pel 


ES col Xs2 (€) 


y(0=1 0] 0 


De forma simplificada se representará el sistema maestro y esclavo, respectivamente, 
del siguiente modo: 


Xm(t) = AmXm(t) + BmnUm(t) 
Ymlt) = CmXmít) 
Xs(t) = Asxs(t) + Bsus(t) 
ys(t) = cgxs(t) 


donde: 
0 1 0 
As Diml B.o= 3 
m 0 EA m == 
Jm m 
5 ol ato = [o 
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0 1 0 
As=lp _*s| B,= z 
A T 
sho Os. 


La estructura del resto de matrices que aparecen en el modelo es la siguiente: 


Km =IKkm1  Km2] 
Ks = [Kks1  Ks2] 
Rm —= [ma Yim2] 
Rs = [Ts1 Ys2] 
G, = 92 


A partir del modelo puede observarse que la señal de control del maestro, u,, (t), y del 
esclavo, us(t), vienen dadas por: 


Um(t) = Keri Xi (E) E RipXs(t > T) 
+ Fm(t) 
4,[£) = Koxst) Sl R5Xp (€ E T) 
+ GaEn(t —T) 


Si se sustituye la señal de control en las ecuaciones de estado (20) y (21) por las 
expresiones (27) y (28) se obtienen las siguientes ecuaciones de estado: 


im (t) = [Am + BmKm]Xm(t) 
+ Bi Rmx (E =T) 
+ Bin Em (0) 

is(t) = [Am + BmKm]Xm(t) 
+ BmRmxXs(t — T) 
+ Bin G2En(t — T) 


Modelado del entorno 


Un NCS en la realización de sus tareas suele interaccionar de algún modo con el 
entorno, por ejemplo agarrando objetos o proyectando algún tipo de material. La 
dinámica de esta interacción debe considerarse en el diseño del sistema de control para 
contemplar las fuerzas de reacción que se generan. Por lo tanto es necesario contemplar 
un modelo del entorno (Buttazzo, 1997). 


bs 
E 
ms 
Ke 
Xs 
e 


Figura 3. Modelo de entorno tipo Kelvin. 


El modelo de entorno más utilizado en robótica es el modelo de entorno tipo Kelvin 
(T.L et al., 2000), Este entorno se representa como una rigidez en paralelo con un 
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elemento de disipación viscoso, ver Figura 3. La ecuación que define su 
comportamiento ante una fuerza externa F es: 


F = k¿Xo(t) + box ¿(t)moXo (t) 


donde k, es la rigidez del entorno, b,¿ es el coeficiente de rozamiento viscoso, m, es la 
masa equivalente, y x¿(t) indica el desplazamiento del entorno. Para incorporar en el 
modelo del NCS mostrado en la Figura 1 la interacción del esclavo con el entorno, será 
necesario que la matriz K, tenga la siguiente estructura: 


ks = [Ks1  Ks2]= [Kks1 =10 $2 Bol 
= [Ks1 $2] = [k, bg] 
=k;—ÍK¿ Ke] 


De este modo la señal que se obtiene al multiplicar el estado del esclavo por la matriz 
K; será (Francis, 2004): 


O 
= Ki(0) — [ko belrs(t) 
= Kixs(t) 
= [Ke bo]| 
= Kixs(t) 

0,(0) 
= Ksxs(t) — (K20,(t) 
+ b2Ó,(0)) 
= ksxs(t) — 4400 


X51 (€) 
X52 (€) 


Por lo tanto, la señal de control de entrada al esclavo será: 


Us(t) = K¿x(E) + RX mlt —T) 
+ G,Fn(t —T) 
= Ksxs(t) — fa(8) 
+ RX l[t — T) + G7En(t 


Si se desea realimentar hacia el maestro la fuerza de contacto del esclavo con el entorno, 
la estructura de la matriz R,,, debe ser la siguiente: 


Kyxs(t) = [mi Tm2] = [KfKe Kfbe] 


donde Kf es la constante de reflexión de esfuerzos considerada. De esta forma la señal 


que se realimenta del esclavo hacia el maestro es: 
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RinXs(t a T) = [Tm1 Tm2] e 2 T) 
0,(t-T) 
=KpKo Kkpb e E 


= kp(k¿0s(t —T) 
+ b,0,(t — T)) 
= ksfa(t —T) 


Como puede comprobarse, con esta estructura de la matriz R,,, se consigue realimentar 
la fuerza de contacto del esclavo con el entorno hacia el maestro. Como se deduce será 
necesario conocer una estimación adecuada de la constante de rigidez k¿ y del 
coeficiente de rozamiento viscoso b¿ para que el modelo del entorno con el que 
interactúa el esclavo reproduzca fielmente el entorno real y, de este modo, el sistema de 
control permita controlar de forma correcta el NCS real (Ferre, 1997). 


Aproximación del retardo 


Como se puede apreciar en las ecuaciones de estado (29) y (30), aparecen términos 
afectados por el retardo existente en la comunicación, en la ecuación (29) aparece el 
estado del esclavo retardado T segundos: x,(t — T'); mientras que en la ecuación (29) 
aparecen el estado del maestro y la fuerza ejercida por el operador retardados T 
segundos: Xm(t — T) y En(t —T). Para eliminar el retardo de estas ecuaciones de 
estado, de forma que no aparezcan términos en función del tiempo retardado T 
segundos, se va a utilizar la aproximación de Taylor de primer orden. El empleo de la 
aproximación de Taylor para eliminar el retardo de las ecuaciones de estado del NCS, y 
poder, de este modo, abordar el diseño del esquema de control utilizando las ecuaciones 
de estado sin la presencia de términos retardados (Ferre, 1997). 


Aproximación de Taylor de una función real f(x): Si f(x) y sus n derivadas son 
continuas en un intervalo que contiene a x = a, se puede aproximar el valor de la 
función en x = a por 


f0=f( +1 a -a) 
+ E (x — ay? 
+ Í e (x — ayi + 


Si se sustituye en la expresión (37) x por t —T y a por t£ se obtiene que el valor de una 
función retardada T segundos puede aproximarse por: 


T? 
FE-D=FO-TFO+P10 


3 
¿E E ES 


De esta forma, las señales retardadas 7 segundos del esquema de control se pueden 


Ciencia e Ingeniería 2015 Vol. (2) N* (2) Velásquez S, Velásquez R. y Romero H. 


MODELADO DE UN NETWORK CONTROL SYSTEM CON RETARDOS EN LA TRANSMISIÓN 


modelar, utilizando la aproximación de Taylor de primer orden, del siguiente modo: 


EE TS) TAE) 
Xt —T) = Xp (E) — Tim (8) 
E (t —T) = Fn(t) — TÉn(t) 


Por lo tanto, si se sustituye las señales retardadas que aparecen en las ecuaciones de 
estado (29) y (30), por las expresiones mostradas en (39) se obtiene: 


Xm(t) = [Am + BmKmIXm(t) + BmRm. 
TB Ri Xt) + BrFn (t) 
Xis(t) = [As + BysK,]x,(t) + BsRsXm(t) 
+TByRsXm(t) mes BsGz Fin (6) 


Si se considera que la fuerza que ejerce el operador es constante, entonces FÉ (t) = 0 
por lo que las ecuaciones (40) y (41) se convierten en: 


Xm(t) = [Am + BmKm]Xm(t) + BmRmA 
BmRmX*s(€) + BmEn(t) 

Xs(t) = [As + BsKs]xs(t) + BsRsXm (8) 
+TB¿RsX mt) => BsG2Fm(t) 


En el siguiente apartado se continuará con el desarrollo de las ecuaciones expuestas en 
este apartado, solo que se orientara su desarrollo a los resultados y su explicación nos 
llevara a la comprensión de las ecuaciones para propiciar una discusión al respecto de 
sus resultados. 


TI. RESULTADOS Y DISCUSIÓN 


En el apartado anterior se presentaron las ecuaciones de estado que se obtenían al 
modelar el retardo utilizando la aproximación de Taylor de primer orden. Las 
ecuaciones que se obtuvieron fueron (42) y (43). En dicho apartado se consideró que la 
fuerza que ejercía el operador era constante, de forma que F,, (t) = 0 ahora bien Si se 
sustituye la expresión (43) en (42) y también se sustituye la expresión (42) en (43) se 
obtienen las siguientes ecuaciones en forma matricial (Yook et al., 2002): 


hsd Y S(A, + B¿K; — TB¿R¿Bm R 
XmCt) Ñ M(BmRm ar TBmRm(As + Bs 
S(BsG> 5 TBsRsBn) 


Flm(B,, —TB,R,,B.G,) 


donde: 
M=(U—-T?BmRmBsRs)” 
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S=(U=T*BR BR) 
De forma simplificada vamos a representar esta ecuación matricial como: 
Xs(t) = Ax(t) + B,¡Fn(t) + BoÉi (0) 


A partir de esta ecuación se va a obtener una ecuación de estado donde no aparezca la 
derivada de la fuerza del operador, sino simplemente aparezca la fuerza del operador, 
así se podrá aplicar al nuevo sistema el método para diseñar el esquema de control 
presentado en la figura 1. Para eliminar la derivada de la fuerza del operador en la 
ecuación (40) se va a realizar el siguiente cambio de variable tal como se indica en 
(Ferre, 1997): 


2(t) = x(t) — BoFin(t) 


Al efectuar las simulaciones considerando tres NCS distintos. El primero corresponde a 
un NCS típico en el que el maestro es menor que el esclavo (NCS I). El segundo 
corresponde a un sistema en el que el maestro y el esclavo son idénticos (NCS II), y por 
último, se ha considerado un NCS en el que el esclavo es menor que el maestro (NCS 
II). Los parámetros del maestro y del esclavo, junto con la constante de reflexión de 
esfuerzos (kf), considerados en cada NCS se muestran en la Tabla 1. 


Tabla 1. Parámetros considerados de cada NCS 


Característica Maestro Esclavo kf 
NES" ze Jn=1.5kgm? Js =2 kg m? 01 
I bm=11 Nm/(rad/s)  bs=15 Nm/(rad/s) ] 
NOS  m=s Jm = 1.5 kg nm? Js =1.5 kg m2 

TI Dm =11Nm/(rad/s)  bs=11 Nm/(rad/s) 

NOS” Jm = 2 kg m? Js =1.5 kg m2 10 
TI bm=15 Nm/(rad/s)  bs=11 Nm/(rad/s) 


Tal como se explicó previamente, para poder diseñar el control de cada NCS, es 
necesario indicar cuál es la dinámica deseada del esclavo y cuál es la dinámica deseada 
del error entre el maestro y el esclavo. Se va a considerar que en todos los casos se 
desea fijar tanto los polos del esclavo como del error en s=-11, es decir, se desea que el 
polinomio característico del esclavo y del error sea (Slotine et al., 1998): 


p(s) = q(s) =s*2+22s + 121 


La razón de ubicar los polos del esclavo y del error en s=-11 se debe a que los 
reguladores que se obtienen con estas dinámicas proporcionan un comportamiento 
satisfactorio de los tres NCS cuando el retardo es de 0.1 s Con los distintos parámetros 
considerados de maestro y esclavo, así como retardo existente y tipos de entornos, se 
pretende validar sobre un amplio abanico de configuraciones, el método de diseño 
aportado en la presente artículo. Para todas las configuraciones se calculará el valor de 
los reguladores según el procedimiento de diseño y posteriormente se simulará el NCS 


con los valores obtenidos. 
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Figura 4. Posición maestro (xm1(1)) y esclavo (xs1(t));. 


En la Figura 4 se muestra la simulación del NCS, considerando que el operador ejerce 
una fuerza constante de 1Nm, cuando el retardo existente en la comunicación es de 0.1 s 
y el entorno con el que interacciona (k. = 10Nm/rad, b¿ = 0.1 Nm/(rad /s)). Como 
puede observarse en las figuras (4) y (5) la posición y velocidad del esclavo 
(51 (€) y Xs2(t)) siguen respectivamente la posición y velocidad del maestro 


Emi (E) Y Xm2 (6). 


tiempo (s 


Figura 5. Velocidad maestro(X,m2(t))y esclavo (xs2(€)) 


En la Figura 6 se muestra la señal de control real del esclavo (a) y la fuerza de 
interacción del esclavo con el entorno (b). Como puede apreciarse en (a), la señal de 
control del esclavo toma el mismo valor final que la fuerza de interacción del esclavo 
con el entorno, de este modo el esclavo no sigue su interacción con el entorno. 
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05 1 15 2 25 3 35 4 45 


Figura 6. La señal de control real del esclavo (a) y la 
fuerza de interacción del esclavo con el entorno (b). 


Si se simula el NCS con los reguladores obtenidos en el caso de que el retardo existente 
en la comunicación es de 0.5 a 0.98 s y el entorno con el que interacciona el esclavo (ke 
= 10 Nm/rad, be = 0.1 Nm/(rad/s)), se obtienen los resultados que se muestran en la 
Figura 7. Como puede observarse, a pesar de la existencia del retardo, la posición y 
velocidad del esclavo siguen la posición y velocidad del maestro. 


— maestro | | 
--- esclavo 


posicion 


a) 


0.005 + 


b) JN 


Figura 7. Simulación del NCS I considerando un retardo de 0.5 a 0.98s y un entorno 
(ke = 10 Nm/rad, be = 0.1 Nm/(rad/s)). (a) Posición maestro (xm1(t)) y Esclavo 
(xs1(t)); (b) Velocidad maestro (x.2(£)) y esclavo (xs,(1)). 


En la Figura 8 se muestra la simulación del NCS, considerando que el operador ejerce 
una fuerza constante de 1 Nm, cuando el retardo existente en la comunicación es de 0.1 
s y el entorno con el que interacciona (ke=10 Nm/rad, be = 0.1 Nm/(rad/s)). Como 
puede observarse en (a) y en (b) la posición y velocidad del esclavo (xs1(t) y xs2(t)) 
siguen respectivamente la posición y velocidad del maestro (xm1(t) y xm2(t)). 
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a) q 


b) pel 


15 2 2.5 3 35 2 45 


Figura 8. Simulación del NCS II considerando un retardo de 0.1 s y un (ke = 10 
Nm/rad, be = 0.1 Nm/(rad/s)). (a) Posición maestro (Xm1(t)) y esclavo (xs1(t)) 
(b) Velocidad maestro (Xm2(t)) y esclavo (xs2(t)). 


Por último, si se simula el NCS con los reguladores obtenidos cuando el retardo es de 
0.5 segundos y el entorno con el que interacciona (ke = 100 Nm/rad, be = 1 Nm/(rad/s)), 
se observa en la Figura 9 que el NCS se vuelve inestable. Es decir, a pesar de diseñar el 
NCS para que el esclavo siga al maestro y, además establecer sus polos en la parte 
izquierda del plano s, no se logra que el sistema sea estable. 


Como se recordará, el NCS diseñado con un retardo de 0.5 segundos, pero con un 
entorno (ke = 10 Nm/rad, be = 0.1 Nm/(rad/s)), no se desestabilizaba, ver Figura 4 y 5. 
Es decir, en el caso de que el entorno de contacto fuese más suave, ha sido posible 
diseñar el esquema de control de forma satisfactoria. 


4 45 


Figura 9. Simulación del NCS II considerando un retardo de 0.5 s y (ke=100 
Nm/rad, be=1 Nm/(rad/s)). (a) Posición maestro (Xm1(t)) y esclavo (xs1(t)), (b) 
Velocidad maestro (Xm(t)) y esclavo (Xs2(t)). 
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En la Figura 10 se muestra la simulación del NCS, considerando que el operador ejerce 
una fuerza constante de 1 Nm, cuando el retardo existente en la comunicación es de 0.1 
s y el entorno con el que interacciona el esclavo es suave (ke = 10 Nm/rad, be = 0.1 
Nm/(rad/s)). Como puede observarse en (a) y en (b), la posición y velocidad del esclavo 
(xs 1(t) y xs2(t)) siguen respectivamente la posición y velocidad del maestro (xm1(t) y 
xm2(t)). 


b) 


£ A £ 1 n a 
0 05 1 15 2 25 3 35 ES 45 


Figura 10. Simulación del NCS [HI considerando un retardo de 0.1 s y (ke = 10 
Nm/rad, be = 0.1 Nm/(rad/s)). (a) Posición maestro (Xm1(t)) y esclavo (xs1(t)), (b) 
Velocidad maestro (xmx(t)) y esclavo (xs2(t)). 


Se va a considerar que ahora el entorno con el que interactúa (ke =100 Nm/rad, be = 1 
Nm/(rad/s)), de forma que el retardo existente en la comunicación sigue siendo de 0.1 
segundos. 


Con estos supuestos, utilizando los reguladores obtenidos según el procedimiento de 
diseño, se obtiene un comportamiento del sistema inestable, tal como se muestra en la 
Figura 11. Al igual que ocurría con el NCS anterior, en el que el maestro y el esclavo 
eran idénticos, al aumentar la rigidez del entorno, el sistema se ha vuelto inestable. En 
este caso el sistema se desestabiliza a pesar de que el retardo existente en la 
comunicación es de 0.1 segundos. Recordemos que el NCS anterior se volvía inestable 
cuando el retardo era de 0.5 segundos. 
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HE —— maestro 
y o. 


E 
> 
posicion 


. f 
b) 3 peta WN 
1 
Figura 11. Simulación del NCS III considerando un retardo de 0.1 s y (ke = 100 


Nm/rad, be = 1 Nm/(rad/s)). (a) Posición maestro (Xmi(t)) y esclavo (xs1(t)), (b) 
Velocidad maestro (xmx(t)) y esclavo (xs2(t)). 


Si se simula el NCS con los reguladores obtenidos en el caso de que el retardo existente 
en la comunicación es de 0.5 s y el entorno con el que interacciona el esclavo es suave 
(ke = 10 Nm/rad, be = 0.1 Nm/(rad/s)), se obtienen los resultados que se muestran en la 
Figura 12. 


Como puede observarse, la posición y la velocidad del maestro y del esclavo convergen, 
pero la convergencia tarda mucho más tiempo en producirse que con el sistema 
diseñado considerando un retardo en la comunicación de 0.1 s y entorno suave, ver 
Figura 11 (ahora el maestro y el esclavo convergen en aproximadamente 3 s, mientras 
que antes convergían en aproximadamente 1.25 s). 


pa] 
Y 
posicion 


de pa pa s de 
o 05 1 15 2 25 3 35 al 45 5 
tiempo (s) 
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b) - | dl A 


002 


Figura 12. Simulación del NCS IM considerando un retardo de 0.5 s y (ke 
= 10 Nm/rad, be = 0.1 Nm/(rad/s)). (a) Posición maestro (Xm1(t)) y esclavo 
(Xs1(t)) (b) Velocidad maestro (xma(t)) y esclavo (xs2(t)). 


IV. CONCLUSIONES 


El método de diseño se ha obtenido a partir de un modelado del NCS en variables de 
estado en el dominio del tiempo en el que se han considerado todas las posibles 
interacciones que pueden existir entre el operador, el maestro, el esclavo y el entorno. 

El método aportado en el presente artículo proporciona un conjunto de ecuaciones cuya 
resolución permite obtener los valores de los reguladores del esquema de control del 
NCS para que el estado del maestro y el del esclavo converjan, consiguiendo de este 
modo que el esclavo siga al maestro. 


Además el método de diseño es capaz de fijar el comportamiento dinámico de la 
convergencia entre el maestro y el esclavo (comportamiento dinámico del error), así 
como el comportamiento dinámico del esclavo. 

Los aportes principales realizadas en este apartado son: 

Modelado de un NCS con retardos en variables de estado contemplando todas las 
posibles interacciones que pueden existir entre el operador, el maestro, el esclavo y el 
entorno. 


Utilización de la aproximación de Taylor de primer orden para obtener un modelo de 
estado del NCS sin la presencia de términos retardados en el tiempo. Método de diseño 
que logra la convergencia de estados entre el maestro y el esclavo, consiguiendo de este 
modo que el esclavo siga al maestro a pesar de la existencia de retardos, y además es 
capaz de fijar tanto la dinámica de la convergencia, como la del esclavo. 

Procedimiento de diseño claro y conciso del sistema de control de NCS con retardos. El 
método de diseño se reduce simplemente a la resolución de un sistema de ecuaciones 
para obtener los parámetros de regulación del esquema de control. 


A partir de las distintas pruebas realizadas, se puede concluir que, en general, los polos 
del polinomio característico del esclavo afectan al comportamiento dinámico del NCS 
del siguiente modo: 


e La posición final del sistema decrece exponencialmente al aumentar la wn del 
esclavo, es decir, al alejar los polos del origen. 
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e La velocidad máxima del maestro y la del esclavo decrecen exponencialmente al 
aumentar la wn del esclavo, o sea conforme los polos del esclavo se alejan del 
origen. 

e El tiempo de establecimiento (ts) de la respuesta del maestro y del esclavo 
disminuyen exponencialmente conforme aumenta la on del esclavo, es decir, a 
medida que se alejan los polos del origen. 

e El tiempo de convergencia entre el maestro y el esclavo disminuye a medida que 
aumenta el valor de wn, es decir, al alejar los polos del origen. 

e El valor absoluto máximo de la señal de control del maestro se incrementa a 
medida que aumenta wn, o sea, al alejar los polos del esclavo del origen. 

e El valor absoluto máximo de la señal de control del esclavo disminuye 
exponencialmente conforme crece «wn, es decir, conforme los polos se sitúan 
más alejados del origen. 


En próximas investigaciones se espera establecer el modelado de sistemas de dos o más 
grados de libertad así como la cuantificación del tiempo de retardo en redes específicas, 
proponiendo un esquema o metodología de control óptima para estas aplicaciones. 
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